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Recently the interest has been revived (see e.g. Refs.[l,2]) to calculate 
of the Debye mass beyond the leading term and to investigate of its gauge 
dependence. This problem is strongly connected with the analogous calcula- 
tions of the infrared n44(0)-limit which are more simple but today they are 
reliably known only up to g'^-order. The next-to-leading term (the one of 
order g^) was found many years ago (at first in paper [3] and then in other 
papers [4,5] and [6]) but till now its accuracy is not confirmed. Unfortunately 
these results are qualitatively different and show that the g'^-term. found for 
the infrared n44-limit is a gauge-dependent quantity and its coefficient (that 
is more essential) is very sensitive to keeping the gauge covariance throughout 
calculations. The problem aggravates when the 1144(0, p) -quantity is calcu- 
lated but namely this limit (as it was shown in Ref. [1,2]) is needed to define 
the Debye screening. However there are many reasons at first to find reliably 
the infrared n44(0)-limit by using the temporal axial gauge since this gauge 
is singled out both for keeping the gauge covariance and for calculating the 
Debye screening. 

The goal of this paper is to derive the nonperturbative equation for the 
infrared n44(0)-limit by operating the standard Green function technique 
within the temporal axial gauge. The obtained equation takes into account 
all perturbative graphs (the ring graphs as well as the vertex corrections) 
and its gauge covariance is guaranteed by exploiting the exact Slavnov- Taylor 
identities to find the nonperturbative three-gluon vertex. This vertex is qual- 
itatively different from the bare one and the derived equation is free from any 
divergencies. The g^-term. is found to be a positive correction to the leading 
one and we compare it with the previously known results. 

It is well-known that the temporal axial gauge is convenient for building 
the nonperturbative schemes since the choice of the gauge vector n^ to be 
parallel to the medium one u^ considerably simplifies the Green function 
technique. The exact polarization tensor (in the axial gauge) is determined 
by only two tensor structures [6] 

Il,Ak) = G U, - ^) + (F - G)B,, , (1) 

and the gluon propagator has a rather simple form 

%(.■) ^ ^ U, - ^) + ^1^^ (2) 



The scalar functions F{k) and G{k) are defined as follows 

Gik) = ^E ^M + 0n44(A:)) , F{k) = ^n44(fc) , (3) 

and they should be calculated through the graph (or another) representation 
for n. Due to a peculiarity of the temporal axial gauge the functions V44{k) 
and T>ij{k) are completely eliminated from the formalism but there is a very 
specific singularity {k\ = 0) which requires a very delicate treatment. The 
exact Slavnov- Taylor identity for the three-gluon vertex function has a rather 
simple form 

r,Tf:^{r, P, q) = ^9r''[V;^\p) - V;^\q)] , (4) 

and namely this fact is a doubtless advantage of the axial gauge. Eq.(4) is 
our main instrument and we exploit it to build the nonperturbative vertex 
function for calculating the infrared n44(0)-limit. We also use Eq.(4) in its 
differential form which allows in many cases (see e.g. Ref.[7]) to define the 
exact infrared limit of the three-gluon vertex in a very convenient manner. 
For example, the infrared r4^'^(— p, 0,p)-limit can be easily found from the 
standard identity 

r4t(-P, 0, p) = -^gr'-^^^f^ , (5) 

which directly results from Eq.(4). This limit being exact has a rather simple 
form 



\t(-p, 0, p) = igr'^ |<5,, [1 + ^] + P,|- 



F{p) 

p2 



Pi , (6) 



and depends on one function which determines the usual representation for 
the P^/(j9)-propagator 

^4-/(p)=-(l + ^)j>,P4. (7) 

Unfortunately the limit (6) is not our case but we shall return back to Eq. (6) 
when the nonperturbative expression for the infrared T^ij {0,p,—p)-liinit is 
discussed. 



The exact graph representation for the gluon polarization tensor is well- 
known (see e.g. papers [6,8]) and contains (in the axial gauge) the standard 
four nonperturbative graphs. However if one considers the n44-components 
only two one-loop nonperturbative graphs are essential since the rest graphs 
(the two very complicated ones) are equal to zero exactly. The analytical 
expression for the first two graphs has a rather simple form and after some 
algebra being performed (by taking into account that F^^^ 
equation for m% (where m|; = 1144(0)) is found to be 



-zgf-f>-T,^^) the 



ml 



9^ 

2/9 



E 

Pi 

E 

P4 



cPp 



2j94 ['Dii{p)Tij4{p, -p, 0)Vji 



(8) 



and we are going to solve it keeping the gauge covariance at each step of 
the calculations. Here all functions (including the vertex one) are exact and 
our main problem is to find the nonperturbative expression for the infrared 
^ijiip, -P, 0)-hmit. 

Unfortunately the exact expression for the infrared Tij^lp, —p, 0)-limit 
(which is not Eq.(6)) lies beyond our possibilities and therefore only its non- 
perturbative ansatz will be represented by following the more general formula 
obtained in Ref.[7]. This formula is found to be 



K-{q,r,p) = -i9rU6,,{r,-p,) 



1 



pz 



'G{r) F{r)r, 



lY* Zi -V^Ai 



_ (G{p) F{p)p\ 
\ p^ p^ p^ 

, / F(r) F{p) 



[(pr)5ij -Pirj](r4-p4) 
\ 1 /F(g) F{r) 



pz 



np^ 7^ 



F{V) \ 



(iif\{q - f) 



q2 \y p. 



|2 

F{r) 






p2 J 



r4P4{r-p)45i 



(9) 



and it is valid for any momentum set including the soft domain. To find 



Eq.(9) the standard inverse gluon propagator is used 



(10) 



and the exact Slavnov- Taylor identities were exploited (see Ref.[7] for details). 
The transversal part of the ri^1{q, r,p)- function is omitted from Eq.(9) since 
it is not essential for what follows. Of course, it is necessary to bear in mind 
to exclude all singularities from Eq. (9) for any momentum going to zero. 

The vertex function (9) easily reproduces the exact formula (6) if one 
momentum goes to zero (at first r^ = and then |r| -^ 0) but it is more 
essential that one can exploit this representation in a more general case to 
find the infrared V^j^iP^ ~P-, 0)-limit. The final result has the form 



rS(i>, -P, 0) = -igr'^ 

'Fip)^ 



+ 2 



pI 



p- 



+ pI 



\ p 

1 d 



IpI "^IpI \ p^ I 




d 



(G{p)\ 
d\v\ \ V' J 



[p^Sij - PiPj] 



(11) 



and we apply this representation for treating Eq.(8). The vertex found is 
qualitatevly different from the bare one and we hope that Eq.(ll) will be 
usefull for many applications. 

All algebra within Eq.(8) is very simple and therefore omitted. The equa- 
tion for m|; is found as follows 



m| 



9^ 
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P4 



d^ 



P 



4pl 



[p2 + G{p)]'' 



(27r)3 
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IpI 



pIi + EM p^ + G{p) 



d (F{py 



9lpl 



p^ 



1 + 



F{p) 
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2pl (dG{p)\ jG{p) F{p)\ 
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[p^ + G{p)\^ 



(12) 



and it is the main subject of the following discussion. Eq.(12) correctly 
summarizes not only the ring graphs but essentially exploits the dressed 
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three-gluon vertex in a nonperturbative manner and it is very probably that 
this equation is exact for the g^-term. Of course, Eq.(12) correctly reproduces 
the leading order for m|; if the functions G(p) and F(p) are omitted 



m. 



9^ 



E 

P4 



d^p d 

(2^)3 dp4 



2P4 
p2 



1 
P4 



(13) 



and the standard regularization for the temporal axial gauge is used 

0. 



11 

Pi:pl 



(14) 



No other terms should be taken into account in Eq.(12) since the rest graphs 
(which usually determine the n^j,-tensor) are exactly equal to zero if the 
n44-quantity is only considered. This fact is due to a simple Lorentz tensor 
structure of the bare r4-vertex function and the specific feature of the tem- 
poral axial gauge where the 1^44-function is eliminated from the formalism. 
The g^-term (as well as all the other ones) should be completely determined 
by Eq.(12) and we found that this equation is free from any divergencies (of 
course if the full sum over p^ is taken into account). 

Our next task is to find the appropriate equation for calculating the g^- 
term on the basis of Eq.(12). This is a pure nonperturbative term and arises 
within Eq.(12) when the soft momentum region is exploited. Due to the 
different infrared behaviour of the functions G(p) and F(p) only the latter 
gives the appropriate contribution to reproduce the ^f^-term and all other 
terms being of (/^-order can be omitted. The final equation for the ^f^-term 
(here the dm^^-term) has a rather simple form 



(5m| 




(15) 



and it can be solved independently from Eq.(12). However, there is a question 
with the first term in Eq.(15) which contains the specific singularity of the 
temporal axial gauge and its analytical behaviour is not clear. Nevertheless 



we insist that this term is equal to zero if only the static n44(0)-limit is used 



ni?(P4 = o,|p|^o) 



3/32 



(16) 



and all calculations are performed in the standard infrared manner (when 
the sum over p4 is replaced by one term with ^4 = 0). Finally the above 
equation is found to be 

^2 



5m\ 



g^N f d^p 1 



(3 J (2^)3 IpI 



d nl^o) 
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nr^o) 



p' 



(17) 



where all functions are known and the integral is calculated in the usual 
manner. Our result has a rather simple form 



m. 



g^N 3 ^g^N 



)3/2j2-2 



(18) 



^3 An' 3 

and it is in an agreement with the results [1,2] and [5] if those are considered 
in the Feyman gauge only. The other results (see Ref. [3,4]) should be checked 
to solve reliably the question of the gauge dependence of the ^r^-term. 

In conclusion we note that the result found for the infrared n44(0)-limit 
is not reproduced through the effective action calculated in a constant back- 
ground field ( see e.g. Refs.[9,10]). This is the case only for the leading term 
but in a general case the correspondence seems to be more complicated. This 
fact needs further investigations to establish finally the status of the ^f^-term 
and its connection with the Debye screening. However there is a hope that 
the latter problem can be solved independently from the scenario with the 
magnetic mass and the question with the magnetic screening is separated (at 
least in the temporal axial gauge) into the outstanding problem. There is 
also a nonperturbative equation [11] 

SN^g"^ f d^p d\ dh 
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P4,q4,r4 
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(27r)3 (27r)3 (27r)3 
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6'-'^\p + q + r) X 



(19) 



but its solution which encounters the infrared divergencies is still uncon- 
firmed. It is important that Eq.(19) is also accessible only in the temporal 
axial gauge and it results from the two rest nonperturbative graphs which 
are equal to zero in the above case. 
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